Roper + Pion decays as a test of the radial resonances 

wavefunct ion 

P. St assart* 

University of Liege, Fundamental Theoretical Physics, 
B5a, Sart Tilman, B-4000 Liege 1, Belgium 

Abstract 

We calculate the decays of baryon resonances to the pion + Roper resonance channel in the frame 
of the flux-tube breaking model. The baryon and meson wavefunctions we use are obtained from a 
linear confinement semi-relativistic constituent quark model. The results, in good agreement with 
the existing data, provide a spectrum independant test on the radial form of the main component 
of the Roper resonance. They shed light on a likely orbital nature for the second delta resonance. 
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I. INTRODUCTION 



The description of the baryon resonances spectrum with the use of a one-gluon exchange 
potential as the source of splitting could not reproduce the ordering of the first two nucleon 
resonances, namely the N(1440) (Roper resonance) and the N(1520) fl]. While Goldstone 
Boson exchange models [2] solve the ordering problem, they cannot account for the complete 
structure of baryon and meson spectra. It has been shown |3| that the Hamiltonian used in 
nonrelativistic or semirelativistic quark models must yield a negative value for the Laplacian 
on at least part of the volume occupied by the resonance. Thus it is a likely possibility 
that models using both OGE and GBE potentials should be used in order to reproduce 
baryon spectra correctly. Spectrum data, however, hardly allow to make a realistic tuning 
of the couplings of these interactions, as most splittings are small and of the same order 
of magnitude as the precision on them. Hence it may be interesting to look for other tests 
on the structure of the states. Decay amplitude calculations can be of most interest from 
this point of view. As the Roper Resonance is the most crucial resonance in studying these 
questions, it is most natural to think of calculating the Roper + pion decays of nonstrange 
baryon resonances. The direct data on the Roper + pion channel look scarse, |4| , as no lower 
bound to the decay widths to that channel is mentionned in the summary table besides the 
one of the radial partner of the roper, the A(1600) resonance. For this reason, that decay 
has not been calculated often. One example of this calculation is that of Capstick and 
Roberts p| who calculate the decay of the higher lying, N>2 band resonances. However, 
using the complete data contained in PDG, we can extract limits on the widths for all three 
and four star baryon resonances, including the N=l,2 bands. It is then worth calculating the 
Roper + pion decays of these states. We do this in the frame of the flux-tube model, using 
wavefunctions for the Roper resonance introduced in p. In the next sections we describe 
the flux-tube model briefly, including the decay mechanism, and the way we treat the broad 
width of the Roper itself. Our results and discussion will be found in the last section. 

II. THE FLUX TUBE MODEL 

The long range attraction between quarks is modelized through the flux tube, which 
corresponds to a linear three body confining potential. The spin-independant part of the 
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Hamiltonian then reads 



H = KE + V(r u r 2 ,r 3 )+E . (1) 

The first term stands for the (semirelativistic) kinetic energy, the last one stands for a 
constant, which can be set from the mass data. The potential consists of a short range 
part, which can be taken as the One Gluon Exchange potential, or as the Goldstone Boson 
Exchange potential, or any combination of them, and of a long range part that reads as the 
sum over all three quarks of the corresponding contributions 

Vt R = «r i4 . (2) 

where r i4 is the distance between the quark i and the point where the flux tubes meet at a 
27r/3 angle. For angles larger than this value that point is identical to the position of the 
corresponding quark. 

In the case of the one gluon exchange potential, the short range part contains a scalar term, 
and spin dependant contributions, namely spin spin, tensor and spin orbit. The short range 
part reads as a sum over the quarks i,j = 1,2,3 with 

V? R = V scal + V ss + V T + V so . (3) 

where V sca i = an d the spin dependent terms read as in ref. For each symmetry 
state, we use the corresponding SU(6) spin flavor wavefunction, while the space part reads 

ipn(p, A) = (j) n (p, A)^ (4) 

where p and A are the Jacobi coordinates and <p n (p, A) ensure the orthogonality of the various 
wavef unctions. For the radial excitation, which is essential in this process, we'll use ((J 

0* = iV(l-a(p 2 + A 2 )t) (5) 

where N is a normalization factor, and a is set so as to insure the orthogonality with the 
ground state. The value k=4 corresponds to the harmonic oscillator case. The hyperfine 
interaction is diagonalized in a truncated space corresponding to the first three levels of 
excitation. The resonances are then to be described by the sum of the products of the 



resulting mixing angles by the corresponding symmetry wavefunctions 
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III. DECAY WIDTHS CALCULATION 



The decay mechanism we use is the breaking of an infinite flux tube, also called the 3 Po 
quark pair creation (QPC) model JtJ, which has been shown to be a good approximation of 
the finite extension QCD-inspired flux-tube breaking 0,0]. The only parameter introduced 
by the model, the breaking amplitude 70, is kept as the value yielding the A(1232) — > 
N + tt width 9]. To ensure we use the correct phase space volume, we use the experimental 
value for the mass of each resonance. This allows us to test wave function and mixing angles 
effects without interference due to specific spectral parameters. Moreover, this is essential in 
obtaining the correct phase space contribution for resonances below the iV(1440)7r threshold, 

hn 

as N(1520). Once the partial wave amplitudes Mj^are obtained |9|,ll0|, they yield the partial 
widths Tig using 

_ 

1 Is — 7777 i T\T7 < J-RJ-irJ-m 1 ^/ J-BJ-ZB > yO) 

TT[2J B + l)M B 

where R,7r,B stand for the Roper resonance, the emitted pion and the decaying baryon, k for 
the momentum, E for the recoil energy, Mb for the baryon mass and the last factor is the 
Clebsch-Gordan coefficient in isospin space. As the Roper resonance itself has a large width 
(350MeV), it is necessary to integrate the weighted value of the partial width over the Roper 
mass interval allowed for the decay process, namely from the Roper production threshold 
M/v + M n up to Mb — M n . The weight a (Mr) corresponds to the Relativistic Breit-Wigner 
distribution 

2T 

a ^ = (M* M W 

?r(r2 (M H -Mno ) } 

where Mr is the nominal value of the Roper resonance mass and T the Roper width. 



IV. RESULTS AND DISCUSSION 



In Table I, we display results for the square root of the width, that is the absolute value of 
the decays amplitudes (in MeV 1 / 2 ). These displayed here have been calculated using k=l in 
the Roper wavefunction. No large variations are observed if one uses values of k up to 4, nor if 
one uses a pure radial symmetry for the Roper resonance. Thirteen resonances are listed, for 
which an upper bound can be extracted from PDG data. Resonances for which there exists 
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TABLE I: Square root of widths (MeV) 1 / 2 



Nominal mass Breit-Wigner PDG upper bound 



N(1520)D13 





1.7 


5.2 


N(1675)D15 


0.4 


2.1 


4.2 


N(1680)F15 


0.2 


1.3 


5.2 


N(1700)D13 


0.4 


1.5 


11.9 


N(1710)P11 


9.5 


12.0 


12.2 


N(1720)P13 


1.5 


1.8 


5.7 


A (162CnS31 


8.6 


8.0 


8.8 


A (1700)D33 


1.8 


2.3 


8.9 


A (1750)P31 


2.1 


2.4 


16.5 


A (1905)F35 


0.7 


1.1 


12.4 


A (1910)P31 


17.6 


14.9 


15.1 


A (1920)P33 


6.9 


6.0 


10.0 


A (1950)F37 


2.4 


3.5 


12.6 



TABLE II: Square root of widths (MeV) 1 / 2 



At nominal mass BW weighted PDG data 

k=l 2.8 11.5 8.5l|5 
k=2 5.9 27.2 
k=4 9.5 35.7 

TABLE III: Square root of widths(Mey) 1 / 2 for ip R =V> (l-a(p 2 + A 2 )) fc/4 

a threshold effect due to the Nr? channel, namely the Sll resonances, are not calculated here 
as the coupled channel calculation is beyond the scope of this article. The comparison with 
data show that all model predictions fall below the experimantal upper bound. Although 
some of these bounds are rather large due to the lack of input data, the result is statistically 
significant. We notice that the effect of the integration with Breit-Wigner weight may be 
quite sizeable for resonance masses up to 1700 MeV, and is not negligible either for higher 
masses, as the width of the Roper is large enough for large variations in the amplitude to 
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take place within the integration interval where the Breit-Wigner weight is not negligible. 

Table II displays the results for the A(1600) resonance, using k=l,2,4 in the radial content 
of both the Roper and A (1600) wavefunctions, with the corresponding mixing angles. One 
can notice that while the k=l case yields a width value that lies within the experimental 
bounds, the k=2 and k=4 cases yield widths values far above the upper bound. It should be 
notice that while the content of the Roper resonance is rather insensitive to the wavefunction 

nn 

used (its radial content ranging from 91 to 96 percent of the total content), |6j,lllj this is not 
the case for its counterpart A(1600). Indeed, while the harmonic limit (k=4) corresponds 
to an 81 percent radial content, the k=l case corresponds to a largely orbital content. We 
checked that a pure radial content yields a very large width, whatever the wavefunction used. 
Hence, within the frame of the flux-tube model, we can affirm that the second A resonance 
should not be considered as a radial state, but rather as a mainly orbital excitation, with 
some coupling (mixing angle of order -0.175 |6||) to the radial symmetry state, while the 
true radial A partner of the Roper resonance is likely to be a higher mass, very large width 
resonance, with a very large coupling to the Roper tt channel. In conclusion, the model 
yields a good agreement with the experimental data for the N=l and N=2 resonances 
decay to Roper tt, for which, however, only large experimental intervals are available. The 
agreement is also good for the A(1600) resonance when we use the form of the spatial radial 
wavefunction (k=l) corresponding to a higher short distance density jf|. This sheds light on 
a likely orbital nature of the A(1600) resonance. This indicates that any hyperfine splitting 
interaction should be compatible with an ordering placing the first radial nucleon resonance 
below the first negative parity orbital nucleon resonances, and placing the first radial delta 
resonance above the first positive parity orbital delta resonances altogether. 
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